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Magnetowave induced plasma wakefield acceleration (MPWA) in a relativistic astro-physical 
outflow has been proposed as a viable mechanism for the acceleration of cosmic particles 
^ . to ultra high energies. Here we present simulation results that demonstrate the viability 

qq ' of this mechanism. We invoke the high frequency and high speed whistler mode for the 

driving pulse. The plasma wakefield so induced validates precisely the theoretical prediction. 
This mechanism is shown capable of accelerating charged particles to ZeV energies in Active 
Galactic Nuclei (AGN). 



1 Introduction 

The origin of ultra high energy cosmic rays (UHECR) is an intriguing question in astrophysics. 
Theories categorized as either "top-down" or "bottom-up" scenarios are proposed to answer 
this question. Each scenario faces its own theoretical and observational challenges x . Since the 
observations from HiRes^and Auger 1 ^ 1 confirm the Greisen-Zatepin-Kuz'min (GZK) suppression 
of the cosmic ray the need for top-down exotic models is reduced. Hence the challenge 

to find a viable "bottom up" mechanism for accelerating UHECR becomes more acute. 

Shocks, unipolar inductors and magnetic flares are the three most potent, observed, "con- 
ventional" accelerators that can be extended to account for ~ ZeV(= 10 21 eV) energy cosmic 
rays^l. Radio jet termination shocks and gamma ray bursts (GRB) have been invoked as sites 
for the shock acceleration, while dormant galactic center black holes and magnetars have been 
proposed as sites for the unipolar inductor acceleration and the flare acceleration, respectively. 
Each of these models, however, presents problems^. Evidently, novel acceleration mechanisms 
that can avoid the difficulties fac ed b y these conventional models should not be overlooked. 

Plasma wakefield accelerators ®^ are known to possess two salient features: (1) The energy 



gain per unit distance does not depend (inversely) on the particle's instantaneous energy; (2) The 
acceleration is linear. These features are essential for a good acceleration efficiency. Although 
high-intensity, ultra-short photon or particle beam pulses that excite the laboratory plasma 
wakefields are not available in the astrophysical setting, large amplitude plasma wakefields can 
instead be excited by the astrophysically abundant plasma "magnetowaves" ^ Protons can be 
accelerated beyond ZeV energy by riding on such wakefields. This attractive concept has never 
been validated through self-consistent computer simulations. In this presentation, we report our 
simulation that confirm this concept^. We also discuss this acceleration mechanism in AGN. 



2 Wakefield Excited by the Magnetowave: Theory and Simulation 

To ensure the linear acceleration, we consider wave modes propagating parallel to the external 
magnetic field. In this case, the eigenmodes are circularly polarized with dispersion relations 
given by^^ 
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where the upper (lower) signs denote the right-hand (left-hand) circularly polarized waves, the 
subscript i denotes the ion species, co p and co c are the electron plasma frequency and the electron 
cyclotron frequency respectively. The right-hand polarized, low frequency solution is called the 
whistler wave which propagates at a phase velocity less than the speed of light. For a sufficiently 
strong magnetic field such that co c 3> co p , the dispersion of the whistler mode becomes more linear 
over a wider range of wavenumbers with phase velocity approaching the speed of light (see Fig.l). 
In this regime the traveling wave pulses can maintain their shape over macroscopic distance, a 
condition desirable for plasma wakefield acceleration. 
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Figure 1: (a) Frequency and (b) phase velocity versus wavenumber for different magnetic fields. The vertical solid 
line is the mean value of the pulse wavenumber chosen for the PIC simulation, and the dashed lines its range. 



The whistler wave packet induces a ponderomotive force along the external magnetic field 
direction^. For convenience, this direction is taken to be along the +z axis. The ponderomotive 
force leads to the plasma wakefield E z which, in the co-moving coordinate £ = z — v g t, reads: 
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where k p = co p /v g and the form factor x(C) is given by 
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with E w (£) the field strength of the whistler wave packet and E\\r its maximum value. Assuming 



the driving pulse frequency is centered around uj and its group velocity v g 
maximum wakefield attainable behind the driving pulse is found to be 
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for a Gaussian driving pulse E W (Q = E w exp(— Q 2 /2a 2 ) with \ = \f^k p a ex.\>(— k 2 a 2 /A)/2, 
ao = eE w j m e cuj the "strength parameter" of the driving pulse, and E w i = m e cuj p /e the "wave- 
breaking" field. We note that Eq. [2] was derived under the non-relativistic approximation for the 
electron motion in the plasma. The relativistic generalization of this equation deserves further 
studies. 

We have conducted computer simulations to study the MPWA process driven by a Gaussian 
driving whistler pulse described above. Our simulation model integrates the relativistic Newton- 
Lorentz equations of motion in the se lf-consistent electric and magnetic fields determined by the 
solution to Maxwell's equations ^ ^ | 13| We used a wavepacket with Gaussian width a = 80A/\/2, 
where A is the cell size taken to be unity, and the wavenumber k = 2-7r/60A. The physical 
parameters uj c /ujp = 6, rrii/m e = 2000 and electron collisionless skin depth c/uj p = 30A were 
taken in the simulations. Other numerical parameters used are: total number of cells in the 
z-direction, L z = 8192A = 273c/uj p , average number of particles per cell was 10, and the time 
step uj p At = 0.1. The fields were normalized by (1/3Q)E W {,. 

We set the maximum amplitude E w = 10, which gives the strength parameter ao = 
eE w /m e cuj = 0.11. The pulse was initialized at zq = 500A = 16.Q6c/uj p . To avoid spuri- 
ous effects, we gradually ramped up the driving pulse amplitude until t = lOOu^ 1 , during which 
the plasma feedback to the driving pulse was ignored. After this time, the driving pulse-plasma 
interaction was tracked self-consistently. As the dispersion relation in this regime is not per- 
fectly linear, there was a gradual spread of the pulse width. Thus x an d E w of the driving 
pulse decrease accordingly. As a result, the maximum wakefield amplitude, E™ x , declined in 
time. Even so, it agrees very well with the theoretical value of I?™ ax ~ 0. 266 ( 1/30) E w b- Fig. 2 
is a snapshot of E x and E z at t = 2?>Quj p l . We note that while the driving pulse continues to 
disperse, the wakefield remains extremely coherent. 
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Figure 2: A snapshot of the plasma wakefield, E z (in black), induced by the whistler pulse, E x (in gray). 



3 MPWA Production of UHECR in AGN 



The MPWA production of UHECR can be effective in AGN. Typically, the jet from an AGN 
extends a very long distance with negligible diverging angle. For an AGN with the central black 



hole mass ~ 10 8 M©, it is reasonable to assume tzagn ~ 10 10 cm -3 and -Bagn ~ 10 4 G in the 
core ll4| jf we f ur ther assume that the AGN luminosity approaches the Eddington limit (~ 10 46 
erg/s) and that the AGN jet size is comparable to the Schwarzschild radius of the central black 
hole, then we find that ao ~ y/lOr] and E w b ~ 10 5 V/cm with r\ the fraction of total energy 
imparted into the magnetowave modes. Since the frequency of magnetowave in this case lies 
in the radio wave region, we assume that the observed AGN radio wave luminosity i s the 
result of total mode-conversion from the magnetoshocks at the same frequency. This then gives 
r\ ~ (10~ 3 — 10~ 4 ) and consequently i?™ ax ~ O(10 2 ) eV/cm from Eq. HI Hence the acceleration 
distance to achieve E ~ O(10 21 ) eV is about 10 pc, which is a tiny fraction of the typical AGN 
jet length. 



4 Summary 

Through PIC simulations, we have confirmed the concept of plasma wakefield excited by a 
magnetowave in the magnetized plasma. We have demonstrated how such a wakefield may 
accelerate particles to ZeV energies in AGN. As a first step, we investigated MPWA in the 
parallel-field configuration. Since both poloidal and toroidal field components are inevitable 
in astro-jets, we will further investigate plasma wakefield excitation and acceleration under 
the cross-field configuration. Besides, we have limited our discussions in the linear regime 
ao -C 1, which is applicable to AGN. However, in other astrophysical settings, such as GRB, the 
magnetowave could be much stronger such that ao 1. We will investigate plasma wakefield 
generations in this regime and therefore explore the MWPA production of UHECR in other 
powerful astrophysical sites. 
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